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Abstract 

Using the methods of loop quantum gravity, we derive a frame- 
work for describing an inflationary, homogeneous universe in a purely 
quantum theory. The classical model is formulated in terms of the 
Ashtekar-Sen connection variables for a general subclass of Bianchi 
class A spacetimes. This formulation then provides a means to de- 
velop a corresponding quantum theory. An exact solution is derived 
for the classical Bianchi type I model and the corresponding semi- 
classical quantum state is found to also be an exact solution to the 
quantum theory. The development of a normalizable quantum state 
from this exact solution is presented and the implications for the nor- 
malizability of the Kodama state discussed. We briefly consider some 
consequences of such a quantum cosmological framework and show 
that the quantum scale factor of the universe has a continuous spec- 
trum in this model. This result may suggest further modifications 
that are required to build a more accurate theory. We close this study 
by suggesting a variety of directions in which to take the results pre- 
sented. 



1 Introduction 

The theory of cosmological inflation pQ has seen continued success since its 
introduction, most notably in its accurate description of the observed fluc- 
tuations of the cosmic microwave background (CMB). Inflation is based on 
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a classical theory of cosmology which is generally thought to be invalid in 
regimes close to the big bang where quantum gravitational dynamics must 
be considered. In this paper, we develop an intrinsically quantum description 
of inflation to try and understand the role of quantum gravity in the early 
universe. Our model consists of a spatially homogeneous universe with scalar 
field matter as this provides the necessary exponential expansion of inflation 
(with an appropriate choice of scalar field potential). 

A theory of quantum cosmology should be based on a more general theory 
of quantum gravity. Such a theory of quantum gravity that has proved very 
successful in recent years is the theory of loop quantum gravity. 1 The start- 
ing point for the general theory of loop quantum gravity is the Ashtekar-Sen 
formulation of general relativity jHl El which greatly simplifies the dynamical 
equations and allows for a well defined quantization of the classical theory. 
While this quantization procedure is a very general and powerful construc- 
tion, we follow a simpler route by imposing homogeneity of the underlying 
spacetime. This allows us to circumvent some of the regularization proce- 
dures that are necessary in the full theory. 

The application of loop quantum gravity to cosmological scenarios has 
already seen many intriguing results in the work of Martin Bojowald and 
others (see jH] for a review of the formulation and the main results). They 
show that the discreteness of the spectrum of the volume operator leads to 
many novel effects including a well-defined evolution through the classical 
big bang singularity In this paper, we use a similar approach to that 
of Bojowald to build our inflationary model with the advantage that we are 
able to find an exact quantum state to describe an inflationary universe. 

The framework for this inflationary quantum cosmological model was 
first introduced in [TUj . In that study, a quantum state for a flat, isotropic 
spacetime was derived and discussed. It was found to be normalizable and 
related to a well-known exact quantum state of the general theory of loop 
quantum gravity, the Kodama state [TT] . The purpose of the following paper 
is to present a generalization of the previous isotropic model to a broader 
class of homogeneous spacetimes and to derive exact quantum states for a 
subclass of these models. 

The starting point for building our model is the classical Hamiltonian 

thorough, technical introductions to loop quantum gravity are available 0E1 as well 
as shorter, less technical reviews [HE], the latter being closely related to the work discussed 
in this paper. 
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framework of ^U] which is summarized here in section |21 In section El we 
specialize this framework to consider only homogeneous universes and discuss 
the procedure of quantization. We then focus on a single class of homoge- 
neous universe, Bianchi type I, in section 0] where a family of exact classical 
and quantum solutions are derived. It is shown in section B~T1 that these solu- 
tions can form a normalizable state and we discuss the possible implications 
for the Kodama state. The specialization of these Bianchi I states to flat, 
isotropic spacetimes is presented in El so as to compare with the analogous 
state of P~0j . We summarize our results in section |H1 and suggest some direc- 
tions for future research. Unless otherwise stated, we work in units in which 
c = h = 1 so that the Planck length is f? p = G. 

2 Classical Hamiltonian Scalar Field Model 

The formalism of loop quantum gravity traditionally begins by describing 
general relativity as a constrained Hamiltonian theory in terms of Ashtekar- 
Sen variables [UJ |7j before applying the quantization procedure of Dirac [12] . 
Since we are interested in modelling inflation by including a scalar field into 
our model, it is possible to recast this classical framework into a more familiar 
Hamiltonian language by using the additional matter degrees of freedom. 
This derivation was first presented in ^Uj and so we briefly summarize the 
procedure here. 

The phase space variables used to describe the geometry of a spatial slice 
of spacetime are a complex, self-dual SU (2) connection A l a and its conjugate 
momentum Ef. Here and throughout we use indices a,b,c,... G {1,2,3} 
as spatial indices and k, . . . G {1, 2, 3} as su{2) indices, the latter being 
raised and lowered by the Kronecker delta 5y. The momentum Ef is the 
densitized spatial triad 

E? = det^Xe- 1 )" (1) 

where e l a is the triad or frame field on the spatial slice and (e _1 )° its inverse. 
The connection can be expressed in terms of the three dimensional spin 
connection of the frame field and the extrinsic curvature K a of the spatial 
slice by 

4 = P a + iejit* (2) 
The canonical relationship between these phase space variables is given by 
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the Poisson bracket 



{4(z), E b j {y)} = iP p tiS i j 6'{x,v) (3) 

where x and y are coordinates on the spatial slice. Since the triad and spin 
connection are both real quantities, it may be necessary to impose the reality 
conditions 

Al + Ai = 2V a (4) 
El- El = (5) 

on the phase space variables where denotes complex conjugation. 

As in conventional inflationary cosmology, we consider gravity coupled 
to a scalar field in a potential V{4>) with conjugate scalar momentum 
7i. The gauge, diffeomorphism, and Hamiltonian constraints can be written 
respectively as 

G t = D a Et = d a E? + e'A^El = 0, (6) 



H a = E%F* b + nd a <f> = 0, (7) 



and 



H = \ 2 ^ k EtE) (F ahk + y^e^) + l -ix 2 + l -E™E } l{d a <P){d b <P) = (8) 

where is the totally antisymmetric, Levi-Civita tensor, F^k is the curva- 
ture of the connection A l a , and the cosmological constant A has been absorbed 
into the potential V(0). 

We choose a decomposition of the spacetime manifold M. into K x S by 
choosing space S to be defined as a constant hypersurface. The limitations 
of such a gauge choice are discussed in [TU] but are not important for the 
results of this paper. It is sufficient to know that such a gauge choice is valid 
during the inflationary phase but may break down after the exponential 
expansion. 

To implement this gauge, we choose a basis of spatial vector fields {d a }^ =1 
such that 

d a (p = (9) 
and demand that this constraint be preserved in time by demanding 

{d a <P, H(N)} = d a (N7r) = (10) 
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where TC(N) = j s N7i is the integrated Hamiltonian constraint that gener- 
ates diffeomorphisms in the time direction. Equation (jlUJI implies that the 
densitized lapse N must take the form 

N= k - (11) 

7T 

where k is a constant, chosen to be k = i~ 2 so that the lapse is dimensionless. 

The diffeomorphism and gauge constraints (J7J) and (jUJ) will be solved 
automatically when we eventually reduce to the homogeneous models and so 
we focus only on satisfying the Hamiltonian constraint (jHJ). In light of the 
time gauge condition (jHJ) we can solve the Hamiltonian constraint for 7r to 
get 



7T 



±A-j 2 e^EtE) (f m + S^^j . (12) 



Combining (|8ji . and (|12j) one can show that the integrated Hamilto- 
nian constraint is 



H ^ = iL^^ P L tk ikm ( Fm + ^ £ ^) ■ (13) 

We use <fi to define a time coordinate 2 

T := fy (14) 

with conjugate momentum 

P:=^/vr (15) 

=> {T, P} = 1 (16) 
so that setting the integrated constraint 7i(N) = in (|13|) implies 

H = P (17) 

where we recognize 



H ■= ±§ I J~^ ijkE ? E * (fM + ^P-e abc E^j (18) 



2 Note that slow-roll inflation 1 is achieved by a decreasing <f> (i.e. 4> is "rolling down 
the hill" of the potential). This means that the forward progression of time corresponds 
to a decreasing T. 
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as the Hamiltonian that generates evolution in the parameter T. Herein we 
will use T instead of <fi and denote the scalar field potential as V(T) to mean 
the original scalar field potential V(4>) evaluated at = T/i^. 

To summarize, by choosing an appropriate coordinate gauge, we have 
been able transform the constrained Hamiltonian theory of general relativity 
into a theory whose dynamics are generated by a traditional Hamiltonian 
functional (|18|) . This will be the starting point for reducing to homogeneous 
spacetimes and, eventually, quantizing to obtain a quantum theory for a 
homogeneous universe. 



3 Bianchi Class A Models 

We now proceed to specialize the Hamiltonian formalism of the previous 
section to consider only the so-called Bianchi class A models. The result will 
be a simpler theory that can immediately be quantized, providing a quantum 
framework to describe a homogeneous, inflationary universe. 

A spatially homogeneous spacetime 3 is characterized by a 3-dimensional 
symmetry group acting transitively on the spatial manifold S. Such a group 
can be described by its algebra with structure constants c I JK which can, in 
general, be written as 

c K IJ = M KL e LI j + 5f I A J] (19) 

where M and Aj are any tensors satisfying 

M IJ Aj = 0. (20) 

The Bianchi class A models are those in which Aj = 0. This class consists 
of 6 models (up to isomorphism) and we further restrict our consideration to 
those 5 whose structure constants can be written in the form 

c'jk = n^e'jx (21) 

where the brackets around an index indicate that it is not to be summed 
over. Herein, we will use the term Bianchi class A to mean this subset of 5 
models. 



3 A complete description of homogeneous spacetimes is given in where the reader 
is referred to for a more detail derivation of the facts stated here. 
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It can be shown jT3| that one can always construct a basis of 1-form fields 
{uj 1 } with dual vector fields {Xj} which are left-invariant under the action 
of the symmetry group (i.e. a basis that does not change when acted upon 
by a group element from the left). Such a basis is uniquely determined by 
demanding that it satisfy the "curl" relations 

duo 1 = —c 1 ~ JK u J A u K (22) 

and that its Lie derivative along timelike vector fields orthogonal to the 
spatial slices vanish. It will prove advantageous to utilize such an invariant 
basis in order to simplify the calculations. 

Following the procedure described in we expand the A and E 

fields in terms of these left-invariant bases. The connection and momentum 
can then be written as 

Ai(x,T) = c (J) (T)A> a / , E?(x,T)= p M(T)AiXf (23) 

where A G SO (3) and the functions {cj} and {p 1 } are constant on a given 
spatial slice and only depend on the time variable (JHJ). Please note that 
is not the most general form of the homogeneous variables but is, in fact, a 
diagonalized form that is well suited to the study of Bianchi class A mod- 
els. The 50(3) matrices are, essentially, a rotation of the invariant basis 
introduced in order to obtain this diagonal form of the connection and mo- 
mentum. These matrices are a statement of the gauge degrees of freedom 
in the A and E variables. For the computations that are to follow, we will 
often use the following properties of SO (3) matrices: 

A/A/ = 5j, eyfcAjAX: = e IJK . (24) 

In order to avoid infinite quantities arising from the (possible) infinitude of 
space, we consider a compact submanifold £ C S of the full spatial manifold 
S. Since the spatial manifold is, by definition, homogeneous, no generality is 
lost by this reduction. We parametrize this submanifold choice by defining a 
dimensionless parameter, R, such that 




where the volume element on S is understood. The £ p 3 normalization is 
simply an arbitrary length scale, chosen to be the Planck length. R can then 
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be considered a parameter of the model that measures the size of our sample 
of space. 

The reduction (|2*3*j) of the full theory to Bianchi class A models reduces 
general relativity from a Hamiltonian theory with an infinite number of de- 
grees of freedom to one with only three degrees of freedom, ci, C2, and C3 
and their conjugate momenta p 1 , p 2 , and p 3 . To determine their canonical 
relationship, we must first define cj and p 1 in terms of the variables A\ and 
E?. A convenient definition is 

01 = l^Xf ( 26 ) 

V 1 = ^3 //^ (/) ^. (27) 

With these definitions we can use the relation Q to compute the Poisson 
bracket 

{ CI (T), ^p\T)} = 5j. (28) 

Hence, we recognize (i? 3 ^ p /z)p 7 (T) as being the true canonical momenta to 
the variables cj(T). Of course, the Poisson bracket between different cj 
variables and between different p J variables respectively vanishes. 

In order to write the diffeomorphism and Hamiltonian constraints in terms 
of the new phase space variables (cj,p J ), we must first calculate the curvature 
of the homogeneous connection 

F = dA+ l -[A ) A}. (29) 

By expanding A = A % Ti in terms of the su(2) basis := — ferj (where Oi are 
the Pauli matrices) and using the identities (}2*2*|) and (}2"^|) we can work out 
the curvature components in the left-invariant basis to be 

F) L = -n^c^^A] + c^c^^jAt (30) 

Here, the constants {n 7 }f =1 are the structure constants of the symmetry 
group defined in (|2*T|). 

We can now use these curvature components and the homogeneous mo- 
mentum field (J23)) to calculate the diffeomorphism constraint (J7J) in the left 
invariant basis. The result is the expression 

Hj = P {J) A{ (-nWc^Af + C( /) c (J) eVA|A fc J ) (31) 
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which, by applying the first relation in (J24|) . can be shown to vanish iden- 
tically for all Bianchi class A models. So, by reducing our field variables 
to an intrinsically homogeneous and diagonalized form, we have also, at the 
same time, automatically satisfied the diffeomorphism constraint. Note that 
this result is not necessarily true for more general forms of the homogeneous 
variables. 

An expression for the Hamiltonian (|18|) in terms of the new phase space 
variables can be obtained in a similar fashion. That is, by substituting the 
curvature (|3T)j) and the components of the homogeneous variables (J2HJ) into the 
Hamiltonian and applying the SO (3) relations ()24|) one obtains the Hamil- 
tonian 



H{ Cl ,p J ;T) = ±2i? 3 



p p c lC2 n 6 1 - + 

V cic 2 3 cic 2 / perm. 

"(32) 

Recall that V(T) is the potential energy function for the scalar field <j> = £ p 2 T, 
kept as an arbitrary function throughout the paper. 

In order to determine the reality condition for the new phase space 
variables, one must, first, expand the three dimensional spin connection in 
terms of the (SO (3) rotated) left invariant basis 

r* = r (J) A}^ (33) 

where Tj can be computed [H)] to be 



Ti — - I n J — - + n K — jt — n 1 ) (no summation on indices) (34) 



)X kP J jP J p K 

2 V'" p J + n W~ n Wf 

where (I,J,K) are even permutations of (1, 2, 3). The resulting reality 
conditions are then 

d + cj = 2Tj (35) 

Pi- Pi = 0. (36) 

Due to the presence of the n 1 parameters in (}3"4*|) . we see that the reality 
conditions are dependent on which Bianchi model is chosen. 

Our classical Hamiltonian theory for a homogeneous universe with a scalar 
field is now completely defined. We have phase space variables ci and p 1 with 
the symplectic structure the reality conditions (jHEJ) and (jSHj) and a time 
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dependent Hamiltonian (|32|) that governs the dynamics. The definitions (|23|) 
then relate the phase space variables back to the Ashtekar-Sen variables 
from which we can then infer all of the information concerning the spacetime 
geometry. Since there are only a finite number of degrees of freedom and no 
constraints needing to be satisfied, one can, in principle, quantize this model 
in a straightforward manner following the traditional algorithm of quantizing 
classical Hamiltonian systems [T7j . 



3.1 Zero-Energy Classical Solution 

Before looking at the quantum theory, it would be very desirable to obtain 
a general solution to this classical Hamiltonian theory. As was seen in the 
calculation of the quantum solution to the flat, isotropic inflationary model 
in [TO], such a classical solution can be crucial in finding a solution to the 
corresponding quantum theory. This will again be the case in the following 
section. 

A first attempt at a general, classical solution may be made via Hamilton- 
Jacobi theory ^H] i n which we seek a Hamilton- Jacobi function S(ci, C2,c^ T) 
such that 

where the constant coefficient in the first relation comes from the non- 
standard Poisson bracket (J2EJ)- Such a function is determined by equa- 
tion (|17|) which now reads as a partial differential equation 



dS 2i 
df ~ L 



— ^c c (n 3 ^ - 1 - tipV(T) 9S ) + Cya 
dc x dc 2 12 \ c i c 2 3i? 3 CiC 2 dc 3 J perm. 



(38) 



where we have chosen the positive root of the Hamiltonian. The determina- 
tion of such a Hamilton- Jacobi function is equivalent to solving the equations 
of motion. 

While we are not able to find a completely general solution to (}38|) for arbi- 
trary values of n 1 , we are able to find a "zero-energy" function Sq(ci, c%, C3; T), 
related to the phase space variables as in (pT7|) . such that H = 0. The de- 
termining differential equation is the vanishing of the right hand side of (|3*Hj) 
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which does indeed vanish if 




where (J, J, K) are cyclic permutations of (1, 2, 3). Such a zero-energy func- 
tion is given as 

<Sb(ci, c 2 , c 3 ; T) = ^cic 2 c 3 - ^n J Cj^ (40) 

plus an irrelevant constant. It is clear that this function does not solve 
the full Hamilton- Jacobi equation (jHEj) since the time derivative of So does 
not vanish for arbitrary potentials V. Of course this would be the general 
solution for the case where V is just the cosmo logical constant with no time 
dependence. 

As we will see in the next section (and as was seen in P-Oj), such zero- 
energy functions are often a useful basis for an ansatz that can be used to find 
a general solution to the Hamilton- Jacobi equation. While such a strategy 
proved useful for the Bianch type I model considered in the next section, 
it is not clear how to successfully proceed with such a strategy for general 
Bianchi class A models. 



4 Bianchi Type I Classical and Quantum So- 
lutions 

We now turn our attention to the simplest, specific class A model: Bianchi I. 
The symmetry group of Bianchi I is the three dimensional translation group 
which is Abelian so the structure constants vanish and we have n 1 = 0. 
Furthermore, because of this commutativity, we may use the coordinate 1- 
forms dx 1 as our left-invariant 1-forms, 

J = dx 1 , (41) 

which clearly satisfy the curl relations (}2*2*|) . 

Since the structure constants vanish for the model under consideration, 
it is clear that the 3 dimensional spin connection (|34j) also vanishes and so 
the reality conditions on the configuration space variables read 

cj = -cj (42) 
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which can trivially be satisfied by defining new variables 



new 
c i 



(43) 



where c new can now be taken to be purely real (i.e. are purely imagi- 
nary in Bianchi I). Herein, we will use, simply, cj to denote the new, purely 
real variables. Of course, given the trivial reality conditions (|3l)j) on p 1 , the 
conjugate momenta can be taken to be purely real. With these new cj vari- 
ables, the Poisson bracket relation ()28|) between the phase space variables 
now reads 



{cj, P J } 



R 3 t 



■si. 



(44) 



We can now adapt the derivation of the previous section to this specific 
model. Using the real variables and n 1 = in the Hamiltonian (|32|) we get 
the Bianchi I Hamiltonian 



H(cj, Pj -T) = 2R 3 



p 1 p 2 c 1 c 2 



gno p 3 

3 cic 2 



+ 



eye. 
perm. 



(45) 



where, again, we have chosen the positive root. Given the slightly different 
Poisson bracket (jH)) for the real variables, a Hamilton- Jacobi function S is 
now defined as 

J_ 1 dS p 



P 



RH V da ' 



OS 
df 



(46) 



and the corresponding Hamilton- Jacobi equation P = H reads 



dS_ 
df 



d^dS_ 
dc\ dc 2 



ClC 2 1 



IpVjT) dS 
3R 3 ciC2 dcs 



eye. 
perm. 



(47) 



To obtain a classical solution for this model, we begin with the zero-energy 
solution (HOD 

3R 3 



and take the ansatz 



So 



3R 3 
LV 



V 



C1C2C3 



cic 2 c 3 (l + u(T)) 



(48) 
(49) 



12 



where u(T) is a complex function of T alone. Upon substituting this ansatz 
into the Hamilton- Jacobi equation (|47|) we find that the c/ variables cancel 
out on both sides and we are left with an ordinary differential equation in u 

V, s 2z, s r- 
u = -(1 + u) + -(1 + u)>/3m (50) 
v ip 

where the dot denotes differentiation with respect to T. Since this is an ordi- 
nary, differential equation, there exists a one (complex) dimensional solution 
space and, therefore, we have obtained a one-parameter family of solutions 
to the classical Bianchi I model via the Hamilton- Jacobi function (|49p . 

Armed with a solution to the classical model we can now proceed to build 
a quantum Bianchi I framework following the usual quantization procedure. 
Since cj and p 1 are real valued, we take the corresponding operators Cj and 
p J to be Hermitian 

(c / )t = c / , (p r )t=^. (51) 

The canonical commutation relations between these operators follow from 
the classical Poisson brackets and are given as 

[cj,p J ] = i{ C/ ,^} = ^/, (52) 

[cj, cj] = i{a, c J } = 0, (53) 
{&, P J ] = i{pi, /} = 0. (54) 

We can represent the vectors in the Hilbert space as functions ^ : M 3 — > C of 
the three configuration space variables c±, C2, and C3. In this representation, 
the phase space operators act as 

c,* = c,*, f to = --L-°-*. (55) 

It is easy to check that these satisfy the above commutation relations. 

The next step in this procedure is to construct a Hamiltonian opera- 
tor corresponding to the classical expression (J45|) . In terms of the classical 
variables, we may order the terms of the Hamiltonian (J45|) as 



H = 2R 3 Cl c 2 c 3 



1 J^a/i £ p V ^\ c y c - 



—p — p 1 p 

c 2 (c 3 ) 2 ci V 3c i c 2 / P erm 



(56) 
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Using the operator representation (|55|) we can construct the corresponding 
operator as 



A 2z 

= -^CiC 2 C 3 



I I <> / 1 ^ g \ cy<-. 



C2(c 3 ) 2 9ci ci <9c 2 V 3i? 3 cic 2 <9c 3 / perm. 



• (57) 



Our goal now is to find a solution ^(ci, c 2 , c 3 ; T) to the resulting Schrodinger 
equation 

i W = H9. (58) 
We introduce the semi-classical state 

H> u := e ^ = e^^i+u) (5Q) 



where S u is the classical Hamilton- Jacobi solution ()49|) . If we denote the 
operator under the square root of (JoTj) as Cl so that 

# = -yC lC2 c 3 VCl (60) 

then it can be shown by direct computation that ^> u is an eigenvector of Cl 
with time dependent eigenvalue: 

tW u = 3w(^— (l + w)J * u (61) 



=*V^* tt = W3u0^(l + u)) *«■ (62) 

Hence, the corresponding action of the Hamiltonian operator ()6()j) can im- 
mediately be determined as 

6R 3 i 

The surprising fact is that the semi-classical state (JoUj) is, in fact, a full 
solution to the Schrodinger equation (|5H|). To see this, we compute the time 
derivative of *f? u to be 

d^ u .3R 3 ( V\ x . , T 

aF = z 4y ClC2C3 ("7 (1 + w)+w l '' (( ' 
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and, upon substitution of the differential equation (|5L)|). this becomes 



:Cic 2 c 3 (l + u)V3u^ u . (65) 



dT £ 2 p V 



Comparison of this with the action of the Hamiltonian operator ()63)1 clearly 
shows that ty u solves Schrodinger's equation (jHEJ)- Hence, equation (fBTJj) gives 
a one parameter family of quantum states for an inflationary universe. 



4.1 Normalizability and the Relationship with the Ko- 
dama State 

A famous solution to the quantum constraint equations of quantum gravity 
with a cosmological constant is the so-called Kodama state ^T] . For a general 
(i.e. non-homogeneous) connection, it can be written as 

V K [A] := ATe- iScs W = jsf e ^h^ Ycs[A] ( 66 ) 
where Yc-s[A] is the Chern-Simons invariant 

Y CS [A] : = ^Tr (a A dA + ^A A A A A J (67) 

and M is a normalization factor depending only on the topology of the mani- 
fold [19] . Similar to the homogeneous, inflationary solution found above, the 
Kodama state is also a semi-classical state corresponding to the Hamilton- 
Jacobi functional 

Scs[A] = I/ cs[A] (68) 

which solves the classical constraint equations of gravity. 

One of the motivations for studying our inflationary model was to try to 
find states similar to the Kodama state (f67)|) in a framework in which the 
state's properties could more easily be calculated. In this section, we inves- 
tigate the similarities between the Kodama state ()66|) and our inflationary 
state (joDj) and discuss the relevant properties of the latter. 

In order to see what the Kodama state looks like in a Bianchi spacetime, 
we need to compute the Chern-Simons invariant (J67j) using the homogeneous 
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connection (J23)) . We do this by explicitly expanding the connection A as AVj 
(recall that T, = — |crj) to obtain 

lfcs = ^(n'c? - 2c 1 c 2 c 3 )u 1 Au 2 Au 3 (69) 

where u 1 are the left-invariant basis 1-forms. Substituting this into the 
Hamilton- Jacobi functional ()68|) we get 

3iR% (\ i 2 \ 
Scs = 4A ( ^ c i - c i c 2 c 3 J (70) 

which is the same as the zero-energy Hamilton- Jacobi function So (I40j) with 
the scalar potential V(T) chosen to be constant. Of course, this similarity 
is not so surprising. The Scs function was chosen so that the Hamiltonian 
constraint vanished while So is defined as that function which causes the 
Hamiltonian (|32|) . which is the integral of the square root of the Hamiltonian 
constraint, to vanish. 

Using this result we can specialize the Kodama state to Bianchi I to 
obtain 

■ c, , .3R 3 l p 

m K = Afe~ lScs = VW^^ C1C2C3 . (71) 

Comparing this with the inflationary state ty u from equation (|59|) we see they 
are the same state up to the time varying function u(T) (with the appropriate 
choice of constant potential V(T)). We also emphasize, once again, that both 
the Kodama state ^ k and our new state ^> u are semi-classical states as well 
as full solutions to their respective theories (the former solves the constraint 
equations while the latter solves Schrodinger's equation). These similarities 
suggest that is a modified Kodama state; it is the analogue of the Kodama 
state when one couples a scalar field to the quantum gravitational field. 

Due to the fact that the configuration space variables ci are real, we may 
take the inner product between two states |\?) and |$) to be 

(tf|$)=/ dc l dc 2 dc 3 ^(ci,c 2 ,c 3 )^(c 1 ,c 2 ,c 3 ) (72) 
Jr 3 

where, as usual, ^(c\, c 2 , C3) := (c\, c 2 , Cs\^) and |ci, c 2 , C3) are the simultane- 
ous eigenvectors of the c/ operators. Using this inner product, it is clear that 
both the Kodama state (fTTj) and the modified Kodama state (joTJj) are both 
continuously normalizable (also called delta function normalizable) because 
they are simply a phase. 
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We now recall the fact there there is a one parameter family of solutions 
u(T) to the differential equation (|5U|). each one of which corresponding to a 
different state If we label each solution u w (T) by its value w := u(T ) 
at some initial time T = To then we may construct (truly) normalizable 
wavefunctions by making "wave packets" of these solutions 

c^T) := J dw f(w)* Uw = J dw f^y^^ 1 ^ (73) 

where f{w) is a normalizable function with support on a suitable interval 
similar to that discussed in PU] 4 . By normalizable, we mean that 

J dw |/H| 2 < oo. (74) 

This type of normalization is not possible in the Kodama state ()71|) as there 
is no such one-parameter family of solutions to take wave packets of. 

Here we see that, by introducing a scalar field into our model, we are 
able to describe the gravitational field (encoded by the configuration space 
variables cj) with respect to the scalar field T w <fi. In this framework, the 
modified Kodama state ([59)1 . which is only continuously normalizable in the 
standard inner product, is able to form a normalizable state by making wave 
packets. 



5 The Spatially Flat, Isotropic Model: A Spe- 
cial Case of Bianchi I 

Before discussing the implications of the above results, we wish to write the 
solution of the previous section specifically for a spatially flat and isotropic 
spacetime. In this case, the real, homogeneous configuration space variables 
are all equal 

ci = c 2 = c 3 =: A (75) 
as are their conjugate momenta 

p 1 = p 2 = p 3 =: E. (76) 

In that paper |10|. it was shown that that the value of the function u(T) is constrained 
such that the metric remain real and maintain a Lorentzian signature. Hence, w is similarly 
constrained to such values. 
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The Hamiltonian (J45)) then reads 



H = 2R 3 A 3 



1 If liV 
3 A> E A> E V-fA> E 



(77) 



where we have written the terms in the order in which they will be written 
in the corresponding Hamiltonian operator. 

We will now write the classical and quantum solutions without much 
commentary as the derivation is completely analogous to that of the previous 
section. The classical Hamilton- Jacobi function is 

S v {A,T) = J ^- ) A\l + v{T)) (78) 
where v is a function of T alone and satisfies the ordinary differential equation 

v = v (l + v) + -(l + v)V3^. (79) 
The representation of the A and E operators can be taken as 

AV = A*, ™ = -~* (80) 
and the corresponding solution to Schrodinger's equation (jSHj) is 

:= e *. = e &W+>) (81) 

which, as in f!59|) . is also the semi-classical state. 

This wavefunction ty v is different from the wavefunction derived in JU] 
using the same framework. The difference is due to the use of different 
operator orderings in the Hamiltonian. The benefit of the solution in this 
paper is that it stems from the more general Bianchi I solution (J59)) whereas 
the solution described in jTU] cannot be generalized. Furthermore, is a 
semi-classical state that is also a full solution to Schrodinger's equation and, 
in this respect, is more closely related to the Kodama state, as was discussed 
in the previous section. The wavefunction of ^U] does not share this property 
but is in fact written as a modified semi-classical state based on the classical 
solution dZHD- 
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5.1 The Spectrum of the Scale Factor 

As a simple example of how one may go about gaining information from 
such a quantum cosmological framework, we can ask what the spectrum of 
eigenvalues is for the operator corresponding to the classical scale factor in 
a flat, isotropic universe. Recalling that E$ is just the densitized triad of the 
spatial manifold, one can easily show that the metric takes the form 

ds 2 = —N 2 dT 2 + EiT^dx 1 ) 2 + (dx 2 ) 2 + (dx 3 ) 2 ) (82) 

where N is the undensitized lapse function described in detail in ^U] . Hence, 
E(T) is the square root of the scale factor. 

Since the operator E is simply a (scaled) derivative operator (J8U|) . its 
eigenvalues are analogous to the familiar momentum eigenvalues from ordi- 
nary quantum mechanics. The eigenvectors of E are given by 

iPe(A) = e lER3£ " A (83) 

and are labelled by the real eigenvalues E. As in ordinary quantum mechan- 
ics, if the domain of A is R (as was assumed in and in this paper), then 
the spectrum of E is continuous. This runs contrary to the prediction of loop 
quantum gravity j2] and, more specifically, to the cosmological predictions of 
Martin Bojowald that the spectrum of volume eigenvalues be quantized. 
If the scale factor had discrete eigenvalues, it would suggest that the universe 
cannot evolve continuously in a classical manner but must proceed by some 
quantum evolution. The fact that we do not find this in our model may be 
a hint that our approach to quantum cosmology is too naive. 

Another possibility that this result may suggest is that the domain of A 
should not be taken as R but as a finite interval. If quantum states were 
required to have periodic boundary conditions on such a finite interval then 
the spectrum would be quantized just as angular momentum is quantized in 
ordinary quantum mechanics. Further study is required in order to resolve 
this issue. 

6 Discussion and Future Projects 

We have shown that the incorporation of a scalar field into the Ashtekar- 
Sen formalism of general relativity allows for a choice of time gauge which 
leads to a standard Hamiltonian theory. By restricting this formalism to a 
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class of homogeneous spacetimes (a subclass of Bianchi class A), we obtain a 
Hamiltonian theory with a finite number of degrees of freedom that describes 
a homogeneous universe filled with a scalar field. The finite number of degrees 
of freedom allows us to straightforwardly quantize to obtain a corresponding 
quantum framework to describe inflationary cosmology 

By looking strictly at the Bianchi I model we are able to find a one- 
parameter family of classical solutions. That is, for every solution u(T) of 
the differential equation ()50|) we obtain a homogeneous spacetime solution 
given by the Hamilton- Jacobi function S u (J49|) . The solutions u(T) also label 
a class of quantum states (jo^j) that are full solutions to the corresponding 
quantum theory. These quantum states are the semi-classical states e tSu 
corresponding to the classical solutions given by S u . We considered flat and 
isotropic spacetimes as a special case of Bianchi I in order to provide an 
alternative solution to that give in [TUj . 

We close our discussion with a brief review of what is left to be done and 
some possible directions this work may go in the future. 

Reality conditions and a physical inner product. In the Bianchi 
class A framework of section El the classical configuration space variables 
Ci were, in general, complex numbers. A quantization of such a complex 
Hamiltonian theory requires a specially constructed inner product so that 
expectation values of real quantities turn out to be real. This problem was 
avoided in the Bianchi I case because the 3-dimensional spin connection van- 
ished and we were able to redefine the theory in terms of purely real variables. 
In such a framework, the standard inner product suffices. For other Bianchi 
models, the spin connection does not vanish identically, the cj cannot be 
taken to be real, and a new inner product will need to be defined. In order 
to have a well-defined quantum framework for the class A models, such an 
inner product will need to be found that is applicable to all models. 

Hermitian ordering of the Hamiltonian. Throughout this paper, 
and in > we have used an operator ordering that leads to a non- Hermitian 
Hamiltonian operator. As is well known, a Hamiltonian must be Hermitian in 
order for probability to be conserved. Hence, it may be required to solve the 
quantum theory using a Hermitian ordering of the Hamiltonian in order to 
obtain a probabilistic interpretation for the wavefunction solution. If this is, 
indeed, required, it is hoped that the techniques used in this paper will lead 
to a similar solution to the Schrodinger equation for the new Hamiltonian. 
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Such a solution has been searched for without any success to date. 

However, given the difficulty in interpreting any wavefunction of the uni- 
verse it is not clear that the ordinary interpretation of quantum mechanics 
can be applied to a theory of quantum cosmology. For this reason, we have 
gone ahead with the non-Hermitian Hamiltonian operator as it leads to an 
exact solution analogous to the Kodama state. Of course, such a solution will 
also be a solution to the Schrddinger equation using a Hermitian Hamiltonian 
at the semi- classical level. 

Numerical analysis. All of the cosmo logical solutions based on the 
framework of this paper have depended on the solution to an ordinary dif- 
ferential equation {e.g. equation (|50p) which cannot be written in closed 
form. Some initial numerical analysis has been done on such equations 10 
but, to fully understand the dynamics of the corresponding quantum solu- 
tions ty u (JSHJ), it will be necessary to better understand the dynamics of u(T). 
It seems that the best way to do this is through computational integration 
of the differential equation. 

Another motivation for the research discussed in this paper was to exam- 
ine the dynamics of a quantum Bianchi IX universe in which the symmetry 
group of space is SU(2). The study of classical spacetimes suggests that 
the dynamics of the spatial geometry of a general manifold close to a sin- 
gularity (such as the big bang singularity) can accurately be modelled by 
a Bianchi IX spacetime [21] and that the evolution of such a Bianchi IX 
model is chaotic (THl 122]- A recent result of the loop quantum cosmology 
program [8] is that loop quantum gravity effects become important close to 
such classical singularities and that such effects curb the chaotic evolution 
of the quantum Bianchi IX model which lead to a more regular evolution of 
the spatial geometry in such regimes [THl I2B] • 

It will be interesting to study the Bianchi IX model in our framework 
to see whether or not the evolution is chaotic close to classical singularities. 
Since an exact wavefunction solution was not found for the Bianchi IX model, 
it may be necessary to evolve an initial state numerically using the Bianchi 
IX Hamiltonian operator (the operator version of equation (|H2*|) with n 1 = 1 
for all J) in order to proceed with this analysis. 

Spatial perturbations and the cosmic microwave background. 

The greatest success of standard, inflationary cosmology is the accurate pre- 
diction of the spectrum of perturbations observed in the cosmic microwave 
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background (CMB) [Tj. Due to the exponential expansion of space during the 
inflationary epoch and the dynamics of quantum scalar fields in such space- 
times, scales well below the Planck length are expanded to macroscopic scales 
in most inflationary models [21] . Since many believe that low energy physics 
no longer applies at such small scales, the precise measurement of the CMB 
may provide the first observational evidence for quantum gravity. Indeed, a 
recent semi-classical analysis j2H| suggests loop quantum gravitational effects 
may leave an observational signature in the CMB. 

The general framework described in this paper may prove to be an advan- 
tageous starting point for a more accurate calculation of the CMB spectrum 
that incorporates the quantum nature of the gravitational field. While very 
little work has been done in this regard, one possible direction is to follow 
an analogous linearization procedure to that described in about one of 
our Bianchi models. In that paper, the authours are able to define an inner 
product which is necessary for any computation. While they find that the 
linearized Kodama state is non-normalizable (in the Lorentzian theory) it is 
hoped that the modifications introduced by the scalar field and the choice 
of time gauge described in this paper will lead to a normalizable state. This 
was indeed the case for the purely homogeneous states, as described in sec- 
tion The successful calculation of the CMB spectrum from a theory of 
quantum gravity would then give us an observable prediction which could 
either support or falsify our approach to quantum cosmology. 
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